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1. Introduction 

Zeta function regularization methods are optimally suited for the calculation of the 
contribution of fluctuations of the vacuum energy, of the quantum fields pervading the 
universe, to the cosmological constant. Order of magnitude calculations of the absolute 
contributions of all fields are known to lead to a value which is off by over hundred and 
twenty orders, as compared with the results obtained from observational fits, what is 
known as the new cc problem. This is difficult to solve and many authors still stick to 
the old problem to try to prove that basically its value is zero with some perturbations 
thereof leading to the (small) observed result (Burgess et al., Padmanabhan, etc.) We 
have also addressed this issue recently in a somewhat similar way, by considering 
the additional contributions to the cosmological constant that may come from the 
possibly non-trivial topology of space and from specific boundary conditions imposed 
on braneworld and other seemingly reasonable models that are being considered in the 
literature (mainly with other purposes too). This kind of Casimir effect would play at 
a cosmological scale. If the ground value of the cc would be indeed zero (and there are 
different hints pointing out towards this), we could then be left with this perturbative 
quantity coming from the topology or boundary conditions and, in particular it could 
be the fact that the computed number is of the right order of magnitude (and has 
the right sign, what is also non-trivial) when compared with the observational value. 
This is proven to be true in some of the aforementioned examples. A further step in 
this approach would be to consider the so-called dynamical Casimir effect or Davies- 
Fulling theory. Although there is no clear understanding of how it should be applied 
in cosmology, some considerations regarding its correct renormalization at laboratory 
scales have been made recently and we will refer to them later. 

The ones above are the physical issues we would like to address ultimately. This 
needs first the heavy mathematics of zeta functions. They will be presented in the first 
part of this work in fair detail. The paper is organized as follows, in correspondence 
with the material presented at the Conference. As a tribute to the actual discoverer 
of the zeta function, namely Leonhard Euler, in this Celebration Year, Sect. 2 recalls 
some essential points that lead him to introduce this function — widely considered to 
be the most important function in Mathematics — with a quick view over the many 
extensions of that concept in the following centuries. In Sect. 3 we describe how the 
concept of zeta function of a pseudodifferential operator has become a decisive tool for 
the regularization of quantum field theories, in special in curved space-time, as clearly 
realized by S. Hawking. This is exemplified in Sect. 4 through the regularization of 
the vacuum fluctuations of a quantum system, under some boundary conditions, with 
a reference to the case of the dynamical Casimir effect (moving boundaries), where 
regularization issues are particularly involved. Finally, Sect. 5 is devoted to the possible 
applications of these results in cosmology, concerning the dark energy issue. 
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2. Euler and the Zeta Function 

There are beautiful accounts on how Euler discovered the zeta function (see, e.g. [HE]). 
The harmonic series 

ff = 1 + I + I + I + I + I + ... (!) 

2 3 4 5 6 K J 

was well known to have an infinite sum. Euler asked himself about the 'prime harmonic 
series' 

11111 
PH = 1 + - + - + - + - + — + ---, 2 
235711' y J 

is it finite or infinite? It is a fact that one cannot split the first series into two, one of 
them being the second, as 

1111 \ (\ 1 1 1 1 \ 
1+ 2 + 3 + 5 + 7 + "j + U + 6 + 8 + 9 + l0 + "j (3) 
and try to show that the second is finite (what would mean the first part is infinite) . So 
Euler considered the function 

C (s = 1 + — + — + — + — + ••• 4 

Provided s is bigger than 1, one can certainly split it up as 

1 1111 \ (\ 1 1 1 1 \ 

Now the idea is to prove that when s approaches 1 the first sum becomes divergent. 
Thus this power s was very useful. 

Making things short, a key step in the whole argument is the celebrated factorization 
of the whole zeta function in terms of prime numbers, namely 

1 1 1 1 1 

f( q\ — y^ y^ y^ y^ y^ . . . ( G ) 

sv ; 1-1/2 S 1 — 1/3 S 1-1/5 S 1-1/7* 1 — 1 /11 s v; 
This comes from the fact that for any prime p and any power s > 1, setting x = 1/p 8 
one has the geometric series 

1 _ 1 1 1 
1 - l/p s ~ + p s+ p 2s + p 3s + 
Euler multiplied together these infinite sums to express his infinite product as a single 
infinite sum as 

klS h q ' 

Pi ■■■Pn n 

with pi, . . . ,p n primes, fcj, . . . , k n positive integers, each such combination occurs exactly 
once and the rhs is just a rearrangement of C( s )- It is widely recognized nowadays that 
Euler's infinite product formula for ((s) marked the beginning of analytic number theory. 

Dirichlet modified the zeta function introduced by Euler. Primes were separated 
into categories, depending on the remainder when divided by k: 

T , \ X(l) , X(2) , X(3) , X(4) , ... 
L(s, x) = — + -y- + -3J- + - 47- + ■ • • > (9) 

where xi n ) is a special function now known as a Dirichlet 'character', that splits the 
primes in the required way. It satisfies the conditions: 
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(i) x{ mn ) = x{ m )x{ n )i f° r an y Tn,n) 

(ii) xip) = x(n + k), Vn; 

(iii) x{ n ) = 0) if n i k have a common factor; 

(iv) = 1- 

Any function L(s, where s is a real number bigger than 1 and \ a character, is known 
as a Dirichlet L-series. The Euler zeta function is the special case with xi n ) = 1 f° r au 
n, another example being xi n ) = M n ) (the Mobius function). 

A very crucial generalization, introduced by Bernhard Riemann, was to allow s 
and x{ n ) to be complex. The celebrated Riemann zeta function, subsequently extended 
by Hurwitz, Lerch, Epstein, Barnes, etc. increased the number and importance of the 
zeta function concept decisively. Many results about prime numbers were proven and 
L-series provide still now a powerful tool for the study of the primes. We should mention 
for completeness that the concept of zeta function has been yet much more extended, 
first to the concept of zeta function of a pseudifferential operator (as we are going to 
see next), but also to the orbits and trajectories in dynamical systems, under the form 
of the Selberg zeta function, the Ruelle, the Lefschetz zeta function, and many others 
that lie outside the scope of this brief summary (Arakelov geometry is one of the most 
active developments right now). In Ref. [2] a directory of all known zeta functions can 
be found (there is even one named after the author of the present article, see also Keith 
Devlin's account there). 

3. The Zeta Function of a Pseudodifferential Operator 

A pseudodifferential operator A of order m on a manifold M n is defined through 
its symbol a(x,£), which is a function belonging to the space S m (M. n x M n ) of C°° 
functions such that for any pair of multi-indexs a, (3 there exists a constant C a ^ so that 
\d£dPa{x,£)\ < C ai 0(l + |f|) m -H The definition of A is given, in the distribution sense, 
by 

Af(x) = (2tt)-" J e*<^>a(x, 0/(0 (10) 

/ a smooth function, / G S, recall S = {f e C°°{R n ); sup x \x^d a f{x)\ < oo, Vet, (3 G M n }, 
S' being the space of tempered distributions and / the Fourier transform of /. When 
a(x, is a polynomial in £ one gets a differential operator. In general, the order m 
can be complex. The symbol of a \I/DO has the form a{x, £) = a m (x,£) + a m _i(x,0 + 
•■• + a m _j(x,0 + ■••> being ak(x, £) = bk{x)^ h . The symbol a(x,£) is said to be 
elliptic if it is invertible for large |£| and if there exists a constant C such that 
K^O -1 ! < C(! + for If I > C- An elliptic ^DO is one with an elliptic symbol. 

Pseudodifferential operators [\l/DO] are useful tools, both in mathematics and in 
physics. They were crucial for the proof of the uniqueness of the Cauchy problem [3] 
and also for the proof of the Atiyah-Singer index formula pi]. In quantum field theory 
they appear in any analytical continuation process (as complex powers of differential 
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operators, like the Laplacian) [5]. And they constitute nowadays the basic starting 
point of any rigorous formulation of quantum field theory j6] through microlocalization, 
a concept that is considered to be the most important step towards the understanding 
of linear partial differential equations since the invention of distributions [TJ. 

3.1. Definition of the Zeta Function 

Let A a positive-definite elliptic \^DO of positive order m G 1, acting on the space of 
smooth sections of E, an n-dimensional vector bundle over M, a closed n-dimensional 
manifold. The zeta function (a is defined as 

{ A (s)=trA- a = Y>7 S , Res>- = s , (11) 

j 

where so = dimM/ordA is called the abscissa of convergence of Ca(s). Under these 
conditions, it can be proven that Ca( s ) has a meromorphic continuation to the whole 
complex plane C (regular at s — 0), provided that the principal symbol of A (that is 
a m (x, £)) admits a spectral cut: Lg = {A G C; Arg A = 9, 9\ < 9 < # 2 } , Spec A n Lg = $ 
(Agmon-Nirenberg condition). The definition of Ca( s ) depends on the position of the 
cut Lg. The only possible singularities of (a(s) are poles at Sk = (n — k)/m, k = 
0, 1, 2, . . . , n — 1, n + 1, . . . . M. Kontsevich and S. Vishik have managed to extend this 
definition to the case when m G C (no spectral cut exists) [Sj. 

3.2. ^DOs on Boundaryless Manifolds 

Let M be a compact n-dim C°° manifold without a boundary, E a smooth Hermitian 
vector bundle over M, A a positive ^DO of positive order m in E, elliptic and selfadjoint 
(admissible). The operator e~ tA , namely e~ tA : f i— > u, is the solution operator for the 
heat equation: d t u + Au = 0, with initial value «| t=0 = /■ 

This operator is traceclass Vt > 0, and as t j it satisfies 

00 oo 

tre"' A ~ J2 a i( A ) t(i ~ n)/m + J2Pk(A)t k logt. (12) 

i=o fc=i 

By Mellin transform: 

1 Z" 00 

Us) = ¥n / e-^t^dt, (13) 

Ca(s) has a meromorphic extension with only possible poles at Sj = (n — j)/m, j G N, 
at most simple at Sj ^ — N, and at most double at Sj G — N. Moreover, 

oij(A) = Res s=s T(s)( A (s), (3 k {A) = Res s= . k (s + k)r(s)( A (s) (14) 

The asymptotic expansion of the heat kernel determines the pole structure of Ca(s), and 
vice versa, (i) If A is a differential operator, then: aij(A) = 0, j odd, = 0, Wk. (ii) 

If A > one still has the same results, but now for A—KerA (subtract DimKer to the 
residue at 0). (iii) If Sj G N, then oij(A) is not locally computable [9| [TO]. 
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3.3. The Zeta Determinant 

Let A a \1/D0 operator with a spectral decomposition: {ipi, where I is some set of 

indices. The definition of determinant p2] starts by trying to make sense of the product 
rTie/^i' which can easily be transformed into a 'sum': lnrjj e /Aj = Sigj hi Aj. From 
the definition of the zeta function of A: Ca( s ) = by taking the derivative at 

s = 0: Ca(0) = — Sig/hiAj, we arrive at the following definition of determinant of A 

mi: 

det c A = exp {-Ca(0)] • (15) 

An older definition (due to Weierstrass) is obtained by subtracting in the series above 
(when it is such) the leading behavior of A; as a function of i, as i —>■ oo, until the series 
^2 ieI In Aj is made to converge. The shortcoming is here — for physical applications — 
that these additional terms turn out to be non-local in general and, thus, they are 
non-admissible in a renormalization procedure |13j . 

In algebraic QFT, in order to write down an action in operator language one needs 
a functional that replaces integration. For the Yang-Mills theory this is the Dixmier 
trace, which is the unique extension of the usual trace to the ideal £( li0 °) of the compact 
operators T such that the partial sums of its spectrum diverge logarithmically as the 
number of terms in the sum: cr N (T) = ^^=0 /J-j = 0(\ogN), fi > \i\ > • ■ ■ The 
definition of the Dixmier trace of T is: Dtr T = liniTv^oo \^rN <T N{T), provided that the 
Cesaro means M(o~)(N) of the sequence in N are convergent as N —>■ oo [remember 
that: M(/)(A) = j^j Ii f( u )^]- Then, the Hardy-Littlewood theorem can be stated 
in a way that connects the Dixmier trace with the residue of the zeta function of the 
operator T _1 at s = 1 (see Connes [Ej): Dtr T = lmi s ^ 1 +(,s — l)( T -i(s). 

3-4- The Wodzicki Residue 

The Wodzicki (or noncommutative) residue [15J is the only extension of the Dixmier 
trace to the \l/DOs which are not in £( 1,0 °). It is the only trace one can define in the 
algebra of ^DOs (up to a multiplicative constant), its definition being: res A = 2 
Res s= o tr(AA~ s ), with A the Laplacian. It satisfies the trace condition: res (AB) = 
res (BA). A very important property is that it can be expressed as an integral (local 
form) res A = J 5 , M tr a,- n (x, £) d£ with S*M C T*M the co-sphere bundle on M (some 
authors put a coefficient in front of the integral: Adler-Manin residue) . 

If dim M = n = — ord A (M compact Riemann, A elliptic, n 6 N) it coincides 
with the Dixmier trace, and one has Res s= iC4(s) = - res A' 1 . The Wodzicki residue 
continues to make sense for \l/DOs of arbitrary order and, even if the symbols aj(x, £), 
j < m, are not invariant under coordinate choice, their integral is, and defines a trace. 
All residua at poles of the zeta function of a \^DO can be easily obtained from the 
Wodzciki residue [TO] . 
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3.5. Singularities of (a 

A complete determination of the meromorphic structure of some zeta functions in the 
complex plane can be also obtained by means of the Dixmier trace and the Wodzicki 
residue. Missing for the full description of the singularities in the above are just the 
residua of all the poles. As for the regular part of the analytic continuation, specific 
methods have to be used (see later). It can be proven that, under the conditions of 
existence of the zeta function of A, given above, and being the symbol a(x,£) of the 
operator A analytic in £ _1 at £ _1 = 0, then it follows that 

Res s=Sfc Us) = - res A~ Sk = - [ tr aZ s n k (x,0 d n ~^. (16) 
m m J S * M 

The proof is rather simple and it can be obtained by invoking the homogeneous 
component of degree — n of the corresponding power of the principal symbol of A, 
obtained by the appropriate derivative of a power of the symbol with respect to at 
= 0, namely 



a_ 



r" (17) 

£-1=0 

Then the proof follows constructively, by easy algebraic manipulation. 
3.6. The Multiplicative Anomaly and its Implications 

Given A, B and AB ^DOs, even if (a, Cb and (ab exist, it turns out that, in general, 
det^(AB) det^A det^-B. The multiplicative (or noncommutative) anomaly (or defect) 
is defined as: 

det c (AB) 



8{A,B) =ln 



&*(«) + ft(0)+&(0). (1* 



_det f A det c 5 

Wodzicki's formula for the multiplicative anomaly [T5J [TTJ, ITS] : 

5(A, B) = ™ { ^ A \ B P A „ a(A, B) := A<***B^. (19) 
v ' ; 2 ordA ordB (ord A + ord 5) ' v ' ; y J 

At the level of Quantum Mechanics (QM), where it was originally introduced by 
Feynman, the path- integral approach is just an alternative formulation of the theory. 
In QFT it is much more than this, being in many occasions the actual formulation of 
QFT [19]. In short, consider the Gaussian functional integration 



(20) 



and assume that the operator matrix has the following structure (being each Ai an 
operator): 




A 

B 



(21) 
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where the last expression is the result of diagonalizing the operator matrix. A question 
now arises. What is the determinant of the operator matrix: det(^4S) or det A ■ det Bl 
This issue has been very much on discussion [20l [21] . 

It is difficult to give a general answer to this question, that is, if it is possible 
to give a universal rule on how to choose the right prescription, and if one can do so 
on mathematical grounds only, without invoking any physical arguments. To start, 
we should not forget that the issue at hand at this level is regularization. This 
means, for one, that there may well be different regularized answers that lead, after 
the corresponding renormalization prescription in each case, to the same renormalized, 
physically meaningful result. But the renormalization process will generically mean 
entering into the physics of the problem in order to choose the right criterion. Thus, 
the answer can in general only be given for the particular example considered. There is 
no space here in order to enter into a more detailed discussion [201 EI]- 

Let us just summarize by pointing out the following. First, that a number of 
serious mistakes and wrong results have appeared in the literature because of forgetting 
about the multiplicative anomaly. Second, that the Wodzicki formula provides a very 
convenient and precise way to calculate the anomaly. Third, that this anomaly turns 
often to be physically meaningful, since it usually (but of course not always) happens 
that the two different regularized results obtained do indeed lead to two different results 
after renormalization [201 EI] (therefore the errors that have been committed in the 
literature, even after going through the whole process of regularization/renormalization 
in a seemingly clean way). Fourth, we know of no mathematically sound prescription 
in order to choose the good regularized answer for the determinant, in general. Maybe 
a better answer to this issue may be given, but it will require further investigation. 

3.7. On Determinants 

Many fundamental calculations of QFT reduce, in essence, to the computation of the 
determinant of some suitable operator: at one-loop order, any such theory reduces in 
fact to a theory of determinants. The operators involved are pseudo differential (\PDO), 
in loose terms 'some analytic functions of differential operators' (such as yl + D or 
log(l + D), but not log/}). This is explained in detail in [22J. It is surprising that 
this seems not to be a main subject of study among mathematicians, in particular 
the determinants that involve in its definition some kind of regularization (related to 
operators that are not trace-class). This piece of calculus falls outside the scope of 
the standard disciplines and even many physically oriented mathematicians know little 
about this. The subject has many things in common with divergent series but lacks any 
reference comparable to the book of Hardy [23J. Actually, this question was already 
addressed by Weierstrass in a way not without problems, since it leads to non-local 
contributions that cannot be given a physical meaning in QFT. For completion, let 
us mention the well established theories of determinants for degenerate operators, for 
trace-class operators in the Hilbert space, Fredholm operators, etc. [24J 
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3.8. The Chowla-Selberg Expansion Formula: Basic Aspects 
From Jacobi's identity for the function 

oo 

9 3 (z, r) := 1 + 2 J2 cos(2nz), q := e l7rT , reC (22) 



n=l 

with 



9 3 (z,r) = -L= e z2 ^6 3 (-\—), (23) 



\f-ir \r r 

or equivalent ly 

OO / 00 2 2 

J2 e ~ {n+z)2t = J j J2 e ~^~ cos(27tnz), z,teC, Re t > 0. (24) 

n=— oo n=0 

In higher dimensions the relevant expression is Poisson's summation formula, profusely 
used by Riemann in his original papers (for recent references see [25J, namely 

being / the Fourier transform of /. An important extension of this theory has consisted 
in the introduction of truncated sums since then neither of these fundamental identities 
is directly applicable [25]. Useful results have been obtained also in these cases, which 
are very important in physical applications, in terms of asymptotic series. 

3.8.1. Extended CS formulas (ECS). Consider the zeta function (with Re s > p/2, A > 
0, Re q > 0) 

r i ~\ ~ s 

CW*)=Z) -{n + cf A{n + c) + q = ^ [Q (« + 3 + <?P (26) 



Hgzp 



new 



where the prime indicates that the point n = is to be excluded from the sum (an 
inescapable condition when c± = ■ ■ ■ = c p = q = 0). We can write 

Q (ft + c) + q = Q(n) + L(n) + q. (27) 

3.8.2. Case q^0{Req>0). Then 

(27r)P/y/ 2 - s r(s - p/2) 2 s / 2+ P/ i+2 ir s q- s / 2+ P/* 



C4,c,g(s 



V / de7Z r(s) VdrtAT(s) 
x ^'cos(27rm • c) (m T A^m) s/2 ~ p/4: Kp/2 _ s (27fy/2qm T A- 1 mj , (28) 

an original expression that we have labeled as [ECSl]. After detailed inspection, it is 
easy to see here that the pole at s — p/2, and its corresponding residue 

Res s=p/2 ( AAq (s) = (det A)' 1 / 2 , (29) 

are explicitly given in the formula, which has in all the following properties. 
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(i) It yields the (analytical continuation of) the multidimensional zeta function in terms 
of an exponentially convergent multiseries, valid in the whole complex plane 

(ii) It exhibits singularities (simple poles) of the meromorphic continuation — with the 
corresponding residua — explicitly. 

(iii) The only condition on the matrix, A, is that it must correspond to a (non negative) 
quadratic form, Q. The vector c is arbitrary, while q is (to start) any non-negative 
constant. 

(iv) K v is the modified Bessel function of the second kind and the subindex in l? x , 2 
means that only half of the vectors m G IT participate in the sum. E.g., if we take 
an index m G IP we must then exclude — m, a simple criterion being: one may 
select those vectors in Z p \{0} whose first non-zero component is positive. 

3.8.3. Case C\ — • ■ • — c p — q — 0. This case is a true extension of CS; we will here 
consider the diagonal subcase only [27] 

,1+s P- 1 



UOO = f^j E (det^.)" 1/2 pV4r* r (s - U2s-j) 

oo , 



(30) 

an expression that truly extends the CS formula and we have labeled as [ECS3d] [27] . 

4. On Zeta Function Regularization 

4-1. Some Considerations on Zeta Regularization 

Regularization and renormalization procedures are essential issues in contemporary 
physics [13J. Among the different methods, zeta function regularization — obtained by 
analytic continuation in the complex plane of the zeta function of the relevant physical 
operator in each case — is one of the most beautiful of all. Use of this procedure yields 
the vacuum energy corresponding to a quantum physical system, with constraints of 
very different nature. The case of moving boundaries seems to present quite severe 
difficulties, though some promising approach to deal with them has appeared [28]. Let 
the Hamiltonian operator, H, of our quantum system to have a spectral decomposition 
of the form (think as simplest case in a quantum harmonic oscillator): {Aj, <y9j}j g /, with 
/ some set of indices (it can be discrete, continuous, mixed, or multiple). The quantum 
vacuum energy is obtained as follows [29J 



Ch/ m (-1), (31) 



E/fi = J2(<Pii (H/fi)<pi) = Tr c F//i = ]T(V/i)- 

where (a is the zeta function corresponding to the operator A, and the equalities are in 
the sense of analytic continuation (since, generically, the Hamiltonian operator will not 
be of the trace class). Actually, this (— trace is no trace in the usual sense. It is highly 
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non- linear, as often explained by the author [30] . Some colleagues are however unaware 
of this fact, which has lead to very serious mistakes and erroneous conclusions in the 
literature. 

The formal sum over the eigenvalues is usually ill defined and the last step involves 
analytic continuation, inherent with the definition of the zeta function itself. Also, an 
unavoidable renormalization parameter, fi, with the dimensions of mass, appears in the 
process, in order to render the eigenvalues of the resulting operator dimensionless, so 
that the corresponding zeta function can actually be defined. For lack of space, we 
shall not discuss those basic details here, which are at the starting point of the whole 
renormalization procedure. The mathematically simple-looking relations above involve 
deep physical concepts, no wonder that understanding them has taken several decades 
in the recent history of quantum field theory. 

4-2. On the Zero Point Energy and the Casimir Force 

In an ordinary QFT, one cannot give a meaning to the absolute value of the zero-point 
energy, and any physically measurable effect comes as an energy difference between two 
situations, such as a quantum field satisfying BCs on some surface as compared with 
the same in its absence, or one in curved space as compared with the same field in flat 
space, etc. This difference is the Casimir energy: Ec = Eq° — Eq = | (tr H BC — tr if). 
But here a problem appears. Imposing mathematical boundary conditions (BCs) on 
physical quantum fields turns out to be a highly non-trivial issue. This was discussed in 
detail in a paper by Deutsch and Candelas [31]. These authors quantized em and scalar 
fields in the region near an arbitrary smooth boundary, and calculated the renormalized 
vacuum expectation value of the stress-energy tensor, to find out that the energy density 
diverges as the boundary is approached. Therefore, regularization and renormalization 
did not seem to cure the problem with infinities in this case and an infinite physical 
energy was obtained if the mathematical BCs were to be fulfilled. However, the authors 
argued that surfaces have non-zero depth, and its value could be taken as a handy 
dimensional cutoff in order to regularize the infinities. Just two years after Deutsch and 
Candelas' work, Kurt Symanzik carried out a rigorous analysis of QFT in the presence of 
boundaries [32J . Prescribing the value of the quantum field on a boundary means using 
the Schrodinger representation, and Symanzik was able to show rigorously that such 
representation exists to all orders in the perturbative expansion. He showed also that 
the field operator being diagonalized in a smooth hypersurface differs from the usual 
renormalized one by a factor that diverges logarithmically when the distance to the 
hypersurface goes to zero. This requires a precise limiting procedure and point splitting 
to be applied. In any case, the issue was proven by him to be perfectly meaningful 
within the domains of renormalized QFT. In this case the BCs and the hypersurfaces 
themselves were treated at a pure mathematical level (zero depth) by using Dirac delta 
functions. 

Not long ago, a new approach to the problem has been postulated [33]. BCs on a 



12 



field, 0, are enforced on a surface, S, by introducing a scalar potential, ex, of Gaussian 
shape living on and near the surface. When the Gaussian becomes a delta function, 
the BCs (Dirichlet here) are enforced: the delta-shaped potential kills all the modes 
of 4> a t the surface. For the rest, the quantum system undergoes a full-fledged QFT 
renormalization, as in the case of Symanzik's approach. The results obtained confirm 
those of pij in the several models studied albeit they do not seem to agree with those of 
[32] . They seem to be also in contradiction with the ones quoted in the usual textbooks 
and review articles dealing with the Casimir effect [33] , where no infinite energy density 
when approaching the Casimir plates has been reported. This has been extended by the 
author using methods of Hadamard regularization, what seems to be a new important 
development in this direction [35J. 

5. Quantum Vacuum Fluctuations, Zeta Regularization, and the 
Cosmological Constant 

5.1. Vacuum Energy Fluctuations and the Cosmological Constant 

The issue of the cc has got renewed thrust from the observational evidence of an 
acceleration in the expansion of our Universe, initially reported by two different groups 
[36J. There was some controversy on the reliability of the results obtained from those 
observations and on its precise interpretation, but after new data was gathered, there 
is now consensus among the community of cosmologists that, in fact, an acceleration 
is there, and that it has the order of magnitude obtained in the above mentioned 
observations [371 EH1 [39]. As a consequence, many theoreticians have urged to try to 
explain this fact, and also to try to reproduce the precise value of the cc coming from 
these observations [30 ] I3T] 132]. 

As crudely stated by Weinberg [33], it is more difficult to explain why the cc is 
so small but non-zero, than to build theoretical models where it exactly vanishes |44|. 
Rigorous calculations performed in quantum field theory on the vacuum energy density, 
pv, corresponding to quantum fluctuations of the fields we observe in nature, lead to 
values that are many orders of magnitude in excess of those allowed by observations of 
the space-time around us. Energy always gravitates [35], therefore the energy density 
of the vacuum, more precisely, the vacuum expectation value of the stress-energy tensor 
{T,xv) = —Sg^v appears on the rhs of Einstein's equations: 

R^u - \g»uR = -8nG(f^ u - £ 9ta ,). (32) 

It affects cosmology: T^ v contains excitations above the vacuum, and is equivalent to a 
cc A = 8nG£. Recent observations yield [36] 

A obs = (2.14 ±0.13 x 1(T 3 eV) 4 ~ 4.32 x HT 9 erg/cm 3 

It is an old idea that the cc gets contributions from zero point fluctuations [37] 

he 

E = — ^2u n , u = k 2 + m 2 /h 2 , k = 2n/A. (33) 
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Evaluating in a box and putting a cut-off at maximum k max corresponding to reliable 
QFT physics (e.g., the Planck energy) 

P ~ H -% k ~ 11)12 W (34) 

Assuming one will be able to prove (in the future) that the ground value of the cc 
is zero (as many suspected until recently), we will be left with this incremental value 
coming from the topology or BCs. This sort of two-step approach to the cc is becoming 
more and more popular recently as a way to try to solve this very difficult issue [48J. 
We have then to see, using different examples, if this value acquires the correct order 
of magnitude — corresponding to the one coming from the observed acceleration in the 
expansion of our universe — under some reasonable conditions. We pursue a quite simple 
and primitive idea, related with the glo bal topology of the universe [19] and in connection 
with the possibility that a faint scalar field pervading the universe could exist. Fields of 
this kind are ubiquitous in inflationary models, quintessence theories, and the like. In 
other words, we do not pretend to solve the old problem of the cc, not even to contribute 
significantly to its understanding, but just to present simple and usual models which 
show that the right order of magnitude of (some contributions to) py which lie in the 
precise range deduced from the astrophysical observations are not difficult to get. In 
different words, we only address here the 'second stage' of what has been termed by 
Weinberg [33] the new cc problem. 



5.2. Vacuum Energy Contribution in Different Models 

5.2.1. Simple model with large and small compactified dimensions. We assume the 
existence of a scalar field extending through the universe and calculate the contribution 
to the cc from the Casimir energy density of this field, for some typical boundary 
conditions. Ultraviolet contributions will be set to zero by some mechanism of a 
fundamental theory. We assume the existence of both large and small dimensions (the 
total number of large spatial coordinates being always three), some of which may be 
compactified, so that the global topology of the universe may play an important role 
[4"9| 1501 I5T| |52| [53] . We know [29] that the range of orders of magnitude of the vacuum 
energy density for common possibilities is not widespread (may only differ by a couple 
of digits) and one can deal with two simple situations: a scalar field with periodic BCs 
or spherically compactified [5_U[55]). The contribution of the vacuum energy of a small- 
mass scalar field, conformally coupled to gravity, and coming from the compactification 
of some small (2 or 3) and some large (1 or 2) dimensions — with compactification radii 
of the order of 10 to 1000 the Planck length in the first case and of the order of the 
present radius of the universe, in the second — lead to values that compare well with 
observational data, in order of magnitude, but with the wrong sign. 



5.2.2. Braneworld models. An important issue in all the previous analysis is the specific 
sign of the resulting force. For scalar fields and the usual compactifications or BCs 
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it is impossible to get the right sign corresponding to the accelerated expansion of 
the universe. However, in braneworld models and others involving supergravitons and 
fermion fields we have been able to prove that the appropriate sign can be obtained 
under quite natural conditions. 

Braneworld theories may hopefully solve both the hierarchy problem and the cc 
problem. The bulk Casimir effect can play an important role in the construction (radion 
stabilization) of braneworlds. We have calculated the bulk Casimir effect (effective 
potential) for conformal and for massive scalar fields [56]. The bulk is a 5-dim AdS or 
dS space, with 2 (or 1) 4-dim dS branes (our universe). The results obtained are quite 
consistent with observational data. A difficulty in this case, however, is the comparison 
of the vacuum energy density obtained in five dimension with the one corresponding 
to four dimensions. Even more, six dimensional models are very on fashion now and 
problems of this kind pop up there too [57] . 

5.2.3. Supergraviton theories We have also computed the effective potential for some 
multi-graviton models with supersymmetry [58]. In one case, the bulk is a flat manifold 
with the torus topology MxT 3 , and it can be shown that the induced cc can be rendered 
positive due to topological contributions [59]. Previously, the case of R 4 had been 
considered. In the multi-graviton model the induced cc can indeed be positive, but 
only if the number of massive gravitons is sufficiently large, what is not easy to fit in a 
natural way. In the supersymmetric case, however, the cc turns out to be positive just 
by imposing anti-periodic BC in the fermionic sector. An essential issue in our model 
is to allow for non-nearest-neighbor couplings. 

For the torus topology we have got the topological contributions to the effective 
potential to have always a fixed sign, which depends on the BC one imposes. They are 
negative for periodic fields, and positive for anti-periodic ones. But topology provides 
then a mechanism which, in a natural way, permits to have a positive cc in the multi- 
supergravity model with anti-periodic fermions. The value of the cc is regulated by 
the corresponding size of the torus. We can most naturally use the minimum number, 
iV = 3, of copies of bosons and fermions, and show that — as in the first, much more 
simple example, but now with the right sign! — within our model the observational values 
for the cc can be approximately matched, by making quite reasonable adjustments of 
the parameters involved. As a byproduct, the results that we have obtained [59] might 
also be relevant in the study of electroweak symmetry breaking in models with similar 
type of couplings, for the deconstruction issue. 
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